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results  of  simulation  using  real  world  demand  data,  which  proved  favorable 
to  this  approach.  It  discusses  a modification  to  current  approaches  for 
determining  cost  of  backorder. 

This  report  was  motivated  by  concern  with  management  of  low  failure 
items.  For  such  items  the  discounted  cash  flow  structure  is  most  appropriate. 


INTRODUCTION 


This  report  is  a sequel  to  an  earlier  paper  by  the  same  author  [5]. 

We  address  a single  product,  continuous  review  model  with  stationary 
Poisson  demand,  fixed  lead  time,  and  backorders.  Such  a model  has  been 
effectively  studied  when  mean  demand  is  known.  However,  we  are  concerned 
with  managing  new  items  for  which  only  a Bayesian  prior  distribution 
on  the  mean  is  available.  As  demand  occurs,  the  prior  is  updated  and 
our  control  parameters  are  revised.  These  include  the  reorder  point 
and  reorder  quantity. 

In  the  earlier  paper,  a heuristic  solution  to  finding  an  optimum 
reorder  point,  given  the  reorder  quantity,  was  presented.  It  was  justified 
analytically,  and  validated  by  simulations  using  randomly  generated  data. 

This  report  introduces  an  alternative  cost  structure,  discounted 
cash  flow,  and  adapts  the  heuristic  approach  to  this  cost  structure 
(Chapter  I) . It  reports  on  the  results  of  simulation  using  real  world 
demand  data,  which  proved  favorable  to  this  approach  (Chapter  II).  It 
discusses  a modification  to  current  approaches  for  determining  cost  of 
backorder  (Chapter  III) . 

This  report  was  motivated  by  concern  with  management  of  low  failure 
items.  For  such  items  the  discounted  cash  flow  structure  is  most 
appropriate. 
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CHAPTER  I 


DERIVATION  OF  COMPUTATIONAL  PROCEDURE 

1.1  Cost  Structure 

Costs  considered  are  purchase  costs  and  backorder  costs.  Each  time 
a buy  of  Q units  is  made,  costs  of  C^  + (u)  (Q)  are  incurred  where  M.  is 
unit  price.  Backorder  costs  are  time  weighted  so  that  each  additional 
day  a unit  is  on  backorder,  a cost  of  C,  is  incurred. 

D 

Costs  are  discounted  to  account  for  the  time  value  of  money  and 
the  probability  of  obsolescence.  The  discount  rate  for  a cost  incurred 
at  time  t is  e In  addition,  looking  at  time  t in  the  future,  there 

is  a probability  of  e ot  the  item  will  not  be  obsolete  by  then.  Hence 
the  discounted,  expected  value  of  a cost  C(t)  which  would  be  incurred  in 
the  absence  of  obsolescence  is  C(t)  e ^ e 0t  = C(t)  e Thus,  impact 

of  time  can  be  accounted  for  by  a discount  factor  a equal  to  (i+o) . 

Holding  costs  are  not  considered  because  the  largest  elements  of  \ 
military  holding  costs  are  financial  and  obsolescence  cost,  and  these  \ 
are  more  correctly  considered  as  discussed  than  as  a linear  holding 
cost  [6].  For  low  failure  items  under  uncertainty,  use  of  linear  holding 
costs  is  particularly  inappropriate.  It  is  quite  possible  for  the  initial  I 
purchase  to  remain  on  hand  for  many  years,  in  which  case,  using  linear  / 
holding  costs,  a total  holding  cost  significantly  greater  than  the  price 
of  the  item  is  computed,  whereas  the  only  cash  outlay  besides  purchase 
is  a small  storage  cost. 

1.2  Basic  Approach 

As  in  the  earlier  paper,  we  find  the  difference  in  cost*  AR,  between 
ordering  at  R,  or  delaying  the  order  until  assets  fall  to  (R-l) . Optimum 
R then  satisfies. 

(1.2.1)  AR*  <_  0 

A(R*+1)  > 0 
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We  think  of  stock  being  issued  in  the  order  in  which  it  is  bought. 
Conceptually,  costs  incurred  in  backordering  a demand  are  associated  with 
the  buy  that  provides  the  unit  which  eventually  satisfies  the  demand,  as 
described  more  fully  in  Silver  [8]  or  Kaplan  [7]. 

Letting 

X - demand  rate 
L - lead  time 

Ba(R,Q|A,L)  - discounted  expected  value  of  backorders  associated 
with  a buy  of  Q units  at  level  R.  Backorders 
are  discounted  back  to  the  time  of  order, 
t - time  until  next  demand 
f(t;X)  - distribution  on  t 

Then 

(1.2.2)  AR(Q,X,L)  - C + (u) (Q)  + C B (R,Q|x,L) 

r k* 

- /°vat  f (t ; X)  [Cp  + (u)  (Q)  + CbBa(R-l,Q|x,L)]dt 

For  example,  if  ordering  at  R means  ordering  at  time  tQ,  then  ordering  at 

R-l  means  ordering  at  time  tQ  + t.  B (R-1,q|x,L)  is  then  the  value  of 

backorders  associated  with  that  order,  discounted  back  to  time  t + t. 

°-at 

For  comparability,  it  must  be  further  discounted  (multiplied  by  e ) 

back  to  time  t . 

o 

Note  that  the  term  in  brackets  under  the  integral  sign  in  (1.2.2) 
is  actually  independent  of  the  variable  of  integration.  Hence 

(1.2.3)  AR(Q,X,L)  = Cp  + (u) (Q)  + CbBa(R,Q| X,L) 

- [C  + (u) (Q)  + C B„(R-1,Q|X,L)][  I"  e~atf(t;X)Jdt 

p Dap 

•k 

To  find  R , we  integrate  (1.2.3)  over  the  prior  on  X.  Computationally, 
this  is  approximated  by  using  a histogram  to  approximate  the  prior: 

/AR(Q,X,L)  dX^rS1  AR(Q,Xi,L)g(Xi) 
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1.3  Integral  Solution 

Since  demand  is  Poisson  distributed,  f(t;A)  is  exponential  with 
mean  1/A,  Hence 


(1.3.1) 


-at  . » _at  -At , A 

f e f(t;A)  = / e Ae  dt  = — — 

o A + a 


1.4  Solution  for  B (R,Q  A,L] 


The  following  is  taken  from  unpublished  work  by  W.  Karl  Kruse. 
Consider  each  of  the  Q units  separately  and  index  them  by  i, 

i ->1,2,...Q.  Let  t be  the  time  until  the  i unit  is  demanded,  measured 

* th 

from  the  time  the  order  is  placed.  If  t^  < L,  then  the  i unit  has  a 

backorder  associated  with  it.  This  backorder  is  cleared  up  at  L.  For 

t^  < L,  the  discounted  time  on  backorder  is: 


(1.4.1) 


t _ -at. 

rL  -at.  i 

/ e dt  = e 

i a 


Now  t^  is  equal  to  the  time  until  the  R+i  demand  since  the  R assets 
existing  at  the  time  of  buy  must  first  be  consumed  before  any  of  the  Q 
units  are  used.  Since  demand  is  Poisson,  t^  has  a Gamma  distribution  with 
parameters  R+i  and  A where  Gamma  is  denoted: 


(1.4.2) 


/ . \ -bx  a-1,  a 

g(x;a,b)  = ex  b 


Using  (1.4.1) 


(1.4.3) 


Ba  (R,Q| A,L) 


Q r “ati  -aL. 

Z f g(t  ;R+i,A)-^ -=£ L dt 

i«=l  ° 1 a 3 


Now,  using  (1.4.2) 


L A 
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-(X+a) t . R+i-1  R+i 

“ati  K 

(1.4.4)  e gCt^jR+i.X)  “ e _i X 

(R+i) 


, R+i 


(X+a) 


R+i 


g(ti;R+i,X+a) 


Substituting  (1.4.4)  into  (1.4.3),  and  denoting  the  cumulative  dis- 
tribution function  of  the  Gamma  by  G(X,a,b), 

(1.4.5)  Bq  (R,Q | X ,L)  = 

Q R+i 

I [fj[+£)R+i  G(L;R+i,X+a)  - e G(L:R+i,X)] 

For  computation  purposes,  we  may  use  the  well  known  equivalence: 


(1.4.6)  G(X;  a,b)  = P(a;bx) 

where  P is  the  Poisson  complementary  cumulative  distribution  function  and 
a is  integer. 


1.5  Histogram  Approximation 

The  prior,  for  convenience,  is  assumed  to  be  Lognormal,  although  in  the 
earlier  paper  we  assumed  Gamma.  We  are  not  being  too  inconsistent  since  for 
most  parameter  values  the  Gamma  and  Lognormal  closely  resemble  each  other. 
For  the  Lognormal  [4], 


(1.5.1) 


f(x;u,o)  = 


ax 


727 


1 2 
exp  ( r—  (log  x - u)  ] 


20 


log 


Var(x>  +1 

2 

E (x) 
x 


u = log  Ex(x)  - a 12 


The  following  additional  results  are  derived  in  the  Appendix. 
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(1.5.2) 


(1.5.3) 


f (x)dx 


g(y)dy 


xf (x) dx 


g(y)dy  • Ex(x) 


where  a' ^ = (log  a^  - u)/o  and  g(y)  denotes  the  Normal  density  for  mean 
0,  variance  1. 

Now,  we  would  like  to  develop  a 10  cell  histogram  to  approximate  the 
Lognormal.  For  given  cell  boundaries  a^,  i **  1 - 11,  we  have  from  (1.5.2) 

£ 

Probability  (i1"^  cell)  = f(x)dx  = G(a^+^)  - G(a^) 

ai 


Value  (ith  cell)  = 


x f(x)dx 


f (x)dx 


G(a’±+1^)  - G(a' ± -°) 

G(a'i+1^  ~ G(ap 


where  G(x)  is  the  cumulative  unit  normal  distribution. 

To  choose  values  for  the  a^,  we  used  the  relationship  a^  = exp(Oa^  + u) , 
i.e.  the  inverse  of  the  relationship  for  going  from  a^  to  a^.  We  always 
chose  for  the  a^,  i = 1 - 11,  these  values: 

-5.0  -3.2  -2.4  -1.6  -.8  0.  +.8  +5.0 


al  a2  a3  a4  a5  a6  a7 all 

Note  that  the  a^  are  actual  cell  boundaries  for  integrating  the  unit  normal. 
The  values  shown  were  found  after  some  experimentation. 
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After  experimenting  with  this  approach  it  was  found  that  the  histogram 

means  were  close  to  the  true  mean,  E (x) . The  histogram  standard  deviations 

th  x th 

(based  on  values  for  i cell,  probability  of  iC  cell,  i = 1 - 11)  were 
dnly  a fraction  of  the  true  standard  deviation.  This  fraction  depended  on 
the  true  coefficient  of  variation, y,  but  given  y,  were  independent  of  the 
true  mean  or  variance.  The  fractional  values  are: 

y .01  0.1  0.2  0.4  0.6  0.8  1.0  1.2  1.4  1.6  1.8  2.0  2.5 

Fraction  .977  .974  .974  .972  .970  .967  .963  .958  .953  .947  .940  .932  .911 

Thus,  if  the  true  y were  1.4,  we  would  find  that  histogram  standard 

deviation  s-  true  standard  deviation  was  .953.  To  partially  correct  for 
this,  if  y were  1.4  we  would  first  divide  the  true  standard  deviation  by 
.953  and  then  calculate  the  histogram  for  this  revised  value. 
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CHAPTER  II 


SIMULATION 

2.1  Purpose  of  Simulation 

The  purpose  of  the  simulation  was  to  determine  whether  the  procedure 
derived  would  work  effectively  in  a real  world  environment,  in  the  manage- 
ment of  a particular  class  of  items.  The  items  of  concern  are  termed 
"insurance."  They  are  highly  essential  with  very  low  failure  rates. 

When  first  introduced  into  the  system,  there  is  no  failure  rate  forecast 
available  for  insurance  items.  Instead,  a prior  distribution  is  developed 
based  on  the  mean  and  variance  of  the  failure  rates  over  the  existing 
catalog  of  insurance  items. 

The  inventory  model  originally  proposed  to  manage  insurance  items  [1] 
was  used  as  the  standard  of  comparison  by  which  to  judge  the  new  procedure. 
This  inventory  model,  in  form,  is  that  appropriate  for  Poisson  demand 
with  known  mean  under  steady  state  conditions  and  linear  holding  and 
backorder  costs,  not  discounted.  Uncertainty  about  the  true  mean  is 
accounted  for  by  using  a Compound  Poisson  distribution  for  demand.  In 
particular,  the  prior  is  assumed  to  be  Gamma,  which  means  that  the  posterior 
distribution  of  demand  about  the  forecasted  mean  is  Negative  Binomial.  This 
model  represented  the  state  of  the  art  when  our  earlier  paper  was  written. 

2.2  Data  Base 

The  data  base  used  was  that  used  to  validate  the  current  policy  [2], 

It  consists  of  four  years  of  history  for  630  insurance  items,  for  the  years 
1967  to  1971,  obtained  from  what  is  now  called  US  Army  Missile  Readiness 
Command . 

There  is  one  important  respect  in  which  the  data  is  treated  differently 
in  this  study  than  it  was  in  the  original  study.  All  requisitions  are 
assumed  to  be  for  one  unit  regardless  of  what  the  historical  quantity  re- 
quested was.  This  "rewriting"  of  history  is  based  on  the  fact  that  insurance 
items  by  definition  should  not  normally  have  requisitions  for  more  than 
one  each.  Currently,  there  is  a recommendation  being  considered  to  establish 
a Maximum  Release  Quantity  of  one  for  insurance  items  except  in  special  cases. 
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Following  are  the  means  and  variances  of  the  catalog  based  on  the 
first  two  years  of  the  data,  computed  with  and  without  the  historical 
requisition  quantity: 


Mean 

Variance 


Historical 

Qtys 

1.85 

28.3  v 


Qtys  Set  To 
1 

.79 

2.1 


The  data  base  has  two  additional  drawbacks  besides  the  problem 
already  described.  Because  of  the  way  it  was  developed,  all  items  in  the 
data  base  were  classified  as  Insurance  at  the  end  of  the  third  year  of 
data;  l.e.,  any  item  which  had  become  obsolete  between  the  first  and  third 
years  would  not  show  up  in  our  catalog  at  all.  Secondly,  for  reasons  not 
known,  the  4th  year  of  data  reflected  a significantly  lower  overall  demand 
rate  as  shown  below: 


Year 

1 

2 

3 

4 

Avg: 
of  3 
Years 

Overall 

Average 

Mean 

o 

00 

• 

.77 

.83 

.39 

.80 

.70 

Variance 

1.84 

2.29 

1.81 

.98 

1.98 

1.73 

2.3  Simulation  Methodology 

All  items  were  started  with  on  hand  equal  to  their  R+Q  values.  The 
Q values  used  were  the  classical  EOQ,  except  that  during  the  first  year  only 
30%  of  EOQ  was  used.  All  of  this  basically  corresponds  to  what  is  currently 
being  done. 

Actual  demands  were  taken  from  the  historical  record,  parameters 
updated  after  each  demand,  and  buys  placed  as  necessary.  Each  item's  pro-r 
curement  lead  time  was  recorded  in  the  data  base  and  these  were  assumed 
to  be  accurate . 

Two  sets  of  mean/variance  pairs  were  used  for  the  prior.  The  3 year 
catalog  averages  were  used  until  an  item  was  2.5  years  old.  Then  the  4th 
year  catalog  averages  were  used  to  reflect  the  drop  off  in  demand  to  be 
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expected  in  the  4th  year.  An  alternative  would  have  been  to  discard  the 
4th  year  entirely,  but  we  felt  it  would  be  more  desirable  to  use  all  the 
data  available. 

It  was  assumed  all  items  had  a unit  price  of  $1000.  The  advantage 
of  this  is  that  no  one  or  two  expensive  items  can  then  dominate  the 
simulation  results.  Previous  work  [1]  had  found  no  correlation  between  unit 
price  and  item  behavior.  If  such  a relationship  did  exist,  but  were  not 
detected,  it  would  undoubtedly  be  that  higher  unit  price  items  had  lower 
demands  than  were  reflected  in  the  catalog.  If  this  were  so,  the  baseline 
policy,  which  tends  to  overstock,  would  be  hurt  most. 

The  cost  structure  assumed  was  that  commonly  used  to  evaluate  alter-  \ 
native  inventory  policies  within  US  Army  Materiel  Development  & Readiness  \ 
Command  [3].  This  structure  is  essentially  consistent  with  the  baseline  j 
policy  but  not  with  the  proposed  policy  and  so  should  favor  the  baseline. 

Costs  include  backorder  costs,  holding  costs  and  "excess"  costs. 
Backorder  costs  are  per  requisition  on  backorder  per  year.  Holding  costs 
are  per  dollar  of  inventory  on  hand  per  year.  Excess  costs  represent  the 
difference  between  the  purchase  price  of  assets  on  hand  or  on  order  at 
the  end  of  the  simulation  and  the  calculated  worth  of  these  assets.  Since 
excess  is  calculated  this  way,  it  is  excluded  from  holding  cost  charges.  \ 
Thus,  MIRCOM's  holding  cost  rate  is  23%  including  10%  for  obsolescence;  ) 
in  charging  holding  cost  in  the  simulation,  a rate  of  13%  is  used.  J 

An  algorithm  is  used  to  calculate  worth  of  end  of  simulation  assets 
for  each  item  [3]: 

a.  Item's  end  of  simulation  demand  rate,  D,  is  estimated.  If  D = 0, 
worth  is  set  to  0. 

b.  End  of  simulation  requirements  objective  (R+Q)  is  computed.  If 
assets  < requirements  objective,  worth  is  set  equal  to  purchase  price 
(no  obsolescence) . 

c.  If  assets  exceed  requirements  objective,  worth  of  all  assets  (A) 
above  requirements  objective  is  calculated  as: 

N i-fl  J"*5  1-fl  N+*5 

(UP)  (D)  z ^ (i^)  + (UP)  (A-ND)  (ij|) 

where  i is  interest  rate,  0 is  obsolescence  rate  and  N is  the  largest  integer 
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t-.5 


. v 1 — 0 

such  that  A >_  ND.  Essentially,  each  item  is  given  a worth  of  (UP)  C^jj/ 
where  t is  the  year  in  which  it  would  be  issued  assuming  constant  demand  rate 
of  D,  no  obsolescence.  Note  that  is  discrete  form  of  e v . In- 

cluded in  0 were  the  small  charges  for  storage  (1%)  and  losses  (2%)  + 10% 
obsolescence.  Since  i ■ 10%,  discount  factor  was  .87/1.1. 

One  problem  in  applying  this  methodology  is  how  to  estimate  D and 
compute  requirements  objectives.  Another  problem  is  that  as  the  data  base 
was  formulated,  there  really  was  little  chance  for  items  to  become  obsolete 
and  no  way  to  distinguish  these  items. 

These  problems  were  solved  in  two  different  ways,  so  we  had  two  dif- 
ferent comparisons  of  results.  In  Method  I,  D was  based  only  on  the  last 
2 years  of  history,  so  if  0,  we  were  effectively  treating  item  as  obsolete 
(this  is  method  for  getting  D used  in  the  past  in  other  contexts  [3]). 

In  Method  II,  D was  based  on  all  four  years  of  history  plus  the  prior. 

No  items  could  be  considered  obsolete  although  an  excess  cost  was  still 
computed  when  assets  exceeded  requirements  objective,  with  obsolescence 
pate  set  to  0;  i.e.  factor  used  was  .97/1.1.  Since  no  obsolescence  rate  . 
was  charged  under  Method  II,  holding  cost  rate  input  to  the  two  inventory 
models,  baseline  and  proposed,  was  only  13%.  Under  both  Method  I and 
Method  II  end  of  simulation  requirements  objective  was  based  on  the  base- 
line model  assuming  Poisson  demand  regardless  of  which  model  was  being 
simulated,  so  that  a given  asset  would  be  costed  out  the  same  way  regardless 
of  inventory  policy. 

Since  backorders  were  not  discounted  for  time  of  occurrence  in  the 
cost  evaluation,  we  did  not  discount  for  interest  in  applying  the  proposed 
procedure.  We  did  discount  for  obsolescence  when  running  under  Method  I. 

Cost  of  backorders  was  assumed  to  be  $23,000  under  Method  I and 
$13,000  under  Method  II.  Under  the  baseline  policies  target  stock  avail- 
ability is  set  to 


1007  - ^old  Cost%) (Unit  Price) 

(Backorder  Cost  Rate) 

as  part  of  the  optimizing  procedure  ■ [3] . Thus,  the  backorder  cost  rates 
assumed  implied  99%  availability  target  under  the  baseline  policy  and  that 
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is  why  they  were  chosen.  For  how  backorder  costs  would  be  computed  in 
the  real  world,  see  Chapter  III. 

2.4  Results 

Shown  below  are  the  cost  comparisons  under  both  Methods  I and  II  for 
baseline  policy  [1]  vs  proposed  policy  (Chapter  I). 


POLICY 

HOLDING  COST 

• EXCESS 

B0  COST 

TOTAL 

Method  I 

Baseline 

1628. 

1758. 

407. 

3793. 

Proposed 

1055. 

1161. 

1488. 

3704. 

Method  II 

Baseline 

1634. 

777. 

232. 

2644. 

Proposed 

1064. 

305. 

826. 

2194. 

Under  Method  I, 

proposed  policy  showed  a 

6%  improvement, 

and  under 

Method  II,  improvement  was  17%.  Because  of  the  data  base  problems,  it  seems 
reasonable  to  conclude  only  that  the  proposed  policy  does  in  fact  work, 
without  trying  to  conclude  just  how  much  better  it  is.  When  runs  were  made 
using  only  one  set  of  catalog  parameters,  based  on  all  4 catalog  years, 
results  were  comparable  under  Method  II,  but  under  Method  I improvement 
narrowed  to  2.4%. 

Just  to  confirm  that  the  difference  between  Method  II  and  Method  I 
was  not  due  to  differences  in  the  backorder  cost  parameter,  Method  II  was 
rerun  with  $23,000  cost  of  backorder.  Results  were  like  those  shown. 

2.5  Use  of  Simple  Poisson 

The  baseline  model  without  any  compensation  for  uncertainty  about 
the  mean  was  tested;  i.e.,  the  simple  Poisson  distribution  was  assumed. 

This  is  easier  to  implement  than  the  learning  model,  and  our  earlier  work 
with  the  learning  model  had  indicated  that  attempting  to  compensate  for 
uncertainty  simply  by  UBe  of  the  compound  Poisson  could  do  more  harm  than 
good. 
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Results  follow: 


Method  I 

Policy 

Holding  Cost 

Excess 

B.O.  Cost 

Total 

Baseline 

1628. 

1758. 

407. 

3793. 

Poisson 

1238. 

1279. 

1120. 

3637. 

Method  II 

Proposed 

1055. 

1161. 

1488. 

3704. 

Baseline 

1634. 

777. 

232. 

2644. 

Poisson 

1244. 

442. 

622. 

2308. 

Proposed 

1064. 

305. 

826. 

2194. 

The  Poisson  fell  between  Baseline  and  Proposed  as  regards  amount  of 
inventory  investment  (note  holding  cost  results) . This  is  partly  because 
the  learning  model  may  react  to  uncertainty  by  lowering  stockage  below 
what  it  would  be  in  the  certain  case,  and  partly  because  the  discounted 
cash  flow  model  will  give  lower  reorder  point  solutions  than  the  steady 
state  model,  under  certainty. 

Poisson  actually  did  better  than  learning  model  in  Method  I,  worse 
under  Method  II.  While  it  is  our  judgement,  taking  into  account  the  results 
and  the  difficulties  with  the  data  base,  that  the, learning  model  is 
superior,  results  are  ambiguous  enough  to  consider  the  simple  Poisson  for 
implementation. 
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CHAPTER  III 


BACKORDER  COSTS 

There  are  two  methods  for  assigning  backorder  costs  used  within 

Dept  of  the  Army,  the  "Lambda"  approach,  and  the  "ERPSL"  approach. 

The  "Lambda"  approach  [9]  is  used  at  wholesale  level.  It  is  assumed 

that  the  backorder  cost  per  requisition  backordered  is  independent  of  the 
* 

item  or  quantity  demanded.  This  cost,  or  Lambda",  is  determined  by  its 
impact  on  projected  budgetary  needs  and  supply  performance.  Higher  "Lambda's" 
results  in  both  higher  dollar  requirements  and  higher  supply  performance. 

A "Lambda"  is  chosen  each  year  which  will  compromise  between  spending  goals 
and  performance  goals. 

The  "ERPSL"  approach  [10]  is  being  implemented  at  retail  for  certain 
low  demand,  highly  essential  items  and  is  being  considered  for  application 
at  wholesale  level  to  insurance  items.  Under  this  approach  the  cost  of 
backordering  a part  is  equated  to  the  cost  of  holding  the  next  higher 
assembly  which  would  be  rendered  useless  if  the  part  were  needed  and  not 
available.  This  type  of  approach  has  been  discussed  in  many  places  in 
the  past  and  has  been  used  by  the  Navy  in  support  of  their  Polaris  sub- 
marine program. 

In  the  remainder  of  this  chapter  we  provide  a rationale  for  the  ERPSL 
approach  within  the  discounted  cost  structure.  We  then  justify  a modifica- 
tion. With  no  loss  of  generality  we  assume  the  next  higher  assembly  is 
the  end  item  as  it  often  is  for  insurance  items. 

* 

Let 

u - marginal  value  per  unit  time  of  having  the  end  item 
operable,  assumed  to  be  constant  over  life  of  item. 

C - cost  of  end  item 

* 

The  Lambda  structure  permits  differentiation  by  items  but  this  is  not 
done  currently  except  that  there  is  a different  Lambda  for  stock  fund 
items  than  for  PA  items. 
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Then,  assuming  there  is  no  production  constraint* 

, 00  —fit- 

(3.1)  C = Qf  ue  'dt 

where  a = ,(0+i) . By  algebra, 

. 2 ) u = C/  e a*dt  » a C 

Thus,  the  value  lost  by  deadlining  the  end  item  for  lack  of  a part  is  at 
a rate  of  a C per  unit  time,  or  unit  price  of  end  item 
holding  cost  rate. 

There  is  a need  to  modify  this  result.  While  the  purchase  price  is 
Incurred  at  time  t = o,  the  use  of  the  item  does  not  begin  until  time  t = 0 
+ PLT,  where  PLT  is  procurement  lead  time.  One  can  argue  about  the  extent 

to  which  the  purchase  cost  is  really  incurred  when  the  end  item  is  ordered 

versus  when  it  is  received,  This  does  not  really  matter  much,  however. 
Since  in  our  reorder  point  models  we  also  assume  the  full  purchase  price  is 
Incurred  when  the  item  is  ordered,  and  the  important  thing  is  that  we  are 
consistent  in  our  assumption. 

We  have  now,  corresponding  to  (3.1) 


00  — (If" 

(3.3)  C = / ue  dt 

PLT 


or,  by  algebra, 


(3.4) 


u = C a e 


a*  PLT 


Using  (3.4)  to  cost  out  a backorder  instead  of  (3.2)  has  the  effect 

a.  PLT 

of  multiplying  the  cost  by  e , For  a = 23%,  PLT  = 1 year,  this 
factor  is  1.26. 
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APPENDIX 


TWO  LOGNORMAL  IDENTITIES 

Notation  used  is  that  of  main  text,  Section  1.5.  Also,  let 
y = log  x where  log  refers  to  natural  logarithm. 

(M)  f*i+1  fOOdx  - exp  [ ~ (log  x - u)2]dx 

ai  ai  X0/2?  2a2 


log  ai+l  r 1 . .2,  y 

/ 1 exp  [ j (y-u)  ] e dy 

loe  al  eyo^T  20 


log  a 


/ 

log  a. 


i+1 


81(y)dy 


where  g^(y)  is  density  for  Normal  random  variable  with  mean  u, 
variance  a2. 


ai+l  log  31+1 

(A2)  ^ x f(x)dx  - log  ^ 


— exp  [ (y-u)2]eydy 

a/ 2tT  2o 


Looking  at  all  exponential  terms: 

[y2  - 2uy  + u2  - 2yo2] 

2o 

= — [(y  - (u+o2)2  - (u+o2)2  + u2)] 

2o 

2 

£y  ” (u+°2)32 

2a 


Substituting  in  (A2) 
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